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In this note, we consider the six-vertex model with domain wall boundary conditions, 
defined on a M x M lattice, in the inhomogeneous case where the partition function 
depends on 2M inhomogeneities Xj and Hk- For a particular choice of the set of Xj we find 
a new determinant representation for the partition function, which allows evaluation of the 
bulk free energy in the thermodynamic limit. This provides a new connection between two 
types of determinant formulae. We also show in a special case that spin correlations on 
the horizontal line going through the center coincide with the ones for periodic boundary 
conditions. 
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1. Introduction 

The six- vertex model was first introduced in Jl[]. It was solved exactly by E. Lieb @ 
and B. Sutherland H in 1967 by means of Bethe Ansatz for periodic boundary conditions 
(PBC). Later the six- vertex model was studied for different boundary conditions @],@,@- 
Domain wall boundary conditions were introduced in 1982 [J7). These boundary conditions 
are interesting because they allow to derive determinant representations for correlation 
functions || and the same boundary conditions help to enumerate alternating sign matrices 
H, IIIJ. Recently the bulk free energy was calculated for these boundary conditions [|TT 



In this paper we show that for special choices of inhomogeneities, one can compute the 
free energy and some correlation functions of the system. This observation might be useful 
because we expect some properties of the model to be independent of the inhomogeneities 
i.e. to depend only on the anisotropy parameter. In the simplest situation, the correlation 
functions coincide with the ones for periodic boundary conditions. 



2. Inhomogeneities and Bethe Ansatz 



In this section we define the inhomogeneous six- vertex model with domain wall bound- 
ary conditions and choose the spectral parameters (inhomogeneities) to satisfy Bethe 
Ansatz equations; this will imply special properties of the partition function. 

We now introduce the notations. Given a M x M square lattice with spectral parame- 
ters Xi and Hk attached to the lines and columns, one defines the usual Boltzmann weights 
a, 6, c to be 

a(A, fx) = sinh(7(A — fx + i/2)) 

b(X,n) = sinh(7(A-/z-z/2)) (2.1) 
c(A, fx) = sinh(z7) 

where 7 is the anisotropy. We have fixed boundary conditions for the external edges: 
horizontal (resp. vertical) external edges are outgoing (resp. incoming). The partition 
function is denoted by Zm({Xj}, {fx k }). Using recursion relations satisfied by the Zm, one 



can show that the following determinant formula holds [T^] , [T^] : 

ni<j, fc <M sinh( 7 (A j - fx k + i/2)) sinh( 7 (A J - fx k - i/2)) 



Z M{{Xj},{fX k }) 



lix<j<j'<M sinh( 7 (A J - X r )) Ui< k <k'<M smh( 7 (^ fc - fi k >)) 

sinh(i7) 



det 

l<j,k<M 



sinh(7(Aj - fi k + i/2)) sinh(7(A J - - fx k - i/2)) 



[2.2) 



In the homogeneous case this represention was used in order to evaluate the bulk free 
energy in the thermodynamic limit [Jll|. In a special inhomogeneous case we shall use 



another determinant represention to evaluate the bulk free energy in the thermodynamic 
limit. Let us define our special case. We choose the spectral parameters in relation to 
the Bethe Ansatz. In order to do that it is convenient to introduce the formalism of the 
Algebraic Bethe Ansatz. The Boltmann weights are encoded into the L-matrix 



L{\) 



/a(A) \ 

6(A) c(A) 

c(A) 6(A) 

V a(A)7 



(2.3) 



We next introduce the monodromy matrix T{\) which is an operator acting on C 2 <8> C 2 
(physical space times auxiliary space) defined by 

T(A; (jli, . . . , fi M ) = L(X - Hm)L{\ - Hm-x) ■ ■ ■ £(A - Hi) (2.4) 

where L(A — Hk) acts on the k th factor of the tensor product in the physical space, and 
the auxiliary space. As an operator on the two-dimensional auxiliary space, T(A) can be 
written as 



C(A) D(A) 

where A, B, C, D are operators on the physical space C 2 . The usual transfer matrix 
corresponding to periodic boundary conditions is defined as 

T(A) = A(A)+D(A) (2.6) 

We shall make use of this operator later. Here, our fixed boundary conditions imply the 
following formal expression for the partition function [JF]: 

^M({Aj},{/ife}) = < | |B(Ai;^i,...,^ m )---B(A m ;/Ui,---,/Wm) | T ) (2.7) 

where | T ) = (o) <8M ( r esp. | I ) = (i) 8 ) is the state with all spins up (resp. down). 

Let us specify inhomogeneities. We choose the {Aj} to be divided into two sets 
{A] = Xj, j = 1, . . . , m 1 } and {A 2 = Aj +m i, j = 1, . . . , m 2 } which each satisfy the Bethe 
Ansatz Equations: 

sinh( 7 (AJ-A^+z)) = il BiDh( 7 (A?-^ + »/2)) 
W sinh( 7 (A« - \$ - z)) 11 sinh( 7 (A« - y, k - z/2)) ( ' ] 
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with a = 1,2. We define the left and right states 

(1| = (t|C(A})...C(Ai l i) 
|2>=B(A?)...B(A^)|t> 



(2.9) 



and the flip operator R to be the operator on the physical space that flips all arrows: in 
terms of the usual Pauli matrices, R = rifcli °%- Now using the overall invariance of 
the monodromy matrix under flip: [T(A),Ra a: ] = (where the additional o~ x acts on the 
auxiliary space), one finds that RB(A)R = C(A) and therefore we can rewrite formula 
( |2.7|) in terms of the states we have defined: 

Z M ({A J },{^}) = <1|R|2> 



(2.10) 

m 2 = M/2. As proven in 
the appendix (A. 3), the Bethe state |2) is an eigenstate of R (with eigenvalue ±1). At 
this point we use orthogonality of Bethe states (A. 2) to conclude that 



We now consider the situation where M is even and m 1 



Zm{{^j}, {M) = ±5 {a 1 },{a 2 }<1|2> 



(2.11) 



The non-zero scalar product (square of the norm) is given by the following formula 
dropping the superscripts: 



'l|l) = (sin 7) 



M/2 



M/2 

n < x m x j 



M / 2 .-,//> x .xx 

nsmh(7(Aj - Aj' + i)) 
, 1 sinh( 7 (A J -A J 0) 



det 

l<j,j'<M/2 



(2.12) 

with the following definitions: a(A) and d(A) are the eigenvalues of A(A) and D(A) acting 
on | t ): 

M 



a(A) = sinh( 7 (A - [i k + i/2)) 
k=i 

M 

d(A) = Y[ sinh( 7 (A -fi k - i/2)) 



(2.13) 



k=i 



and the tpj are the logarithms of the B.A.E. ( |2.8| ): 

M/2 



<Pj = i lo g( a ( A i)/ d ( A i)) + i Yl log ~ 



sinh(7(Aj — Aj' + i)) 
sinh(7(Aj — Aj/ — i)) 



[2.U) 



Note that the general determinant formula (|2.2|) , in the case of 2 identical sets {A)} = {A 2 }, 
becomes 

Y\i<j<m/2 sinh 2 (7(Aj - (i k + i/2)) sinh 2 (7(Aj - [x k - i/2)) 



Zm({^}, {^k}) 



Yll<j<j'<M/2 sinh (l( X i - V)) rii<fc<fc'<M sinh( 7 (^ fc - (i k ,)) 

T<k<M 

with </>(A) = sinh(ry)/(sinh(7(A + i/2)) sinh(7(A - i/2))) and ^(A) = ^0(A). Our new 
determinant representation ( |2.12|) is quite different from Eq. ( |2.15|) ; it is in particular much 
easier to study its thermodynamic limit. 



3. Correlation functions 

We have studied so far the function partition of the system; what about correlation 
functions? 

It is unfortunately not possible to find such a simple expression for an arbitrary 
correlation function; however, if we restrict ourselves to the case of special correlations 
which lie on a fixed horizontal line, then one can again derive determinant formulae for 
them. 

To be specific, let us assume again that the {Aj} consist of two identical sets {Xj, j = 
1 . . . M/2}. Then the probability that all arrows located at columns hi, . . . , k n and between 
lines M/2 and M/2 + 1 are up, is given by 



(Tfci •••T*J = , , | z— r — r— tt , =— , , ^ x L (3.1 



Ul 


B(A, 


i) . . . B(A M /2) "i, .••"'*„ B(Ai) . . .B(A A/ / 2 ) 




Ul 


B(A 1 )...B(A M/2 )B(A 1 )...B(A AJ/2 )| 


T> 



where txu = §(1 + = (o o) ac ^ s on ^ e ^ th s P ace - 

Similarly to what was done for the partition function, one can transform Eq. (|3? 
using the flip operator R and find: 

(flfci ■ ..tj = - ■ ' " ( 3 - 2 ) 

In other words this is simply the usual correlation functions of the spin operators ^(l + a z ) 
for the corresponding spin chain. The computation of these averages is a well-known 



problem, and the general strategy to perform it is by now well understood [EHl . One must 
use the solution of the quantum inverse scattering problem for these operators [|15| , |T6|| : 



k-i 



M 



]"]T(// / + //2)A(// A . + //2) ]"[ T(m+i/2) 



(3.3) 



1=1 



l=k+l 



and then use the fact that our state 1 1) is an eigenstate of the T(A), as well as the commu- 
tation relations satisfied by the A(A) and B(A). In this way we can reduce Eq. (|3.1| ) (resp. 
Eq. (|3.2|) ) to a sum of expressions of the type ( | | B(Ai) . . . B(Am/2)B(Ai) . . . B(Am/2) | T ) 
(resp. ( | | C(Ai) . . . C(Am/2)B(Ai) . . . B(Aa//2) | T ))• These can finally be expressed as 
determinants either using Eq. fl2.2| ), or according to the general formula for scalar products 
(p. 237 of §). 

The case n = 1 is trivial: |l) is an eigenstate of R, and therefore (iTk) = 1/2. The 
simplest higher correlation function is the emptiness formation probability where all ki are 
nearest neighbors, in which case we obtain immediately 



(■K k +1 ■ ■ ■ TTfc+n) 



n 

l=k+l 



M/2 



^< + */ 2 ) n s 



inh( 7 (/i/ - Xj + i/2)) 



fj{ sinh(7(^ ; - Xj - i/2)) 
We shall give the thermodynamic limit of this expression in a particular case 



fc+n 

ii n A( fH +i/2)\i) 

l=k+l 

(3.4) 



4. Thermodynamic limit 



We shall show in a particular example how to take the thermodynamic limit in for- 
mulae ( |2.12| ) and (|3T4T). We set all fik to and consider the critical regime i.e. 7 real. 
We specify the state |l) to be the ground state of the transfer matrices T(A) (or of the 
Hamiltonian of the corresponding XX Z spin chain). In the limit M — >• 00 the Xj form a 
continuous distribution on the real axis determined by its density p(A) = 1/(2 cosh(7rA)). 

Let us first consider the free energy. One can show that the undetermined sign in 
Eq. ( 2.11 ) is +, and therefore we have F = — log (l|l)- We now analyze Eq. ( |2 . 1 2 ) in the 
large M limit. We see that we have 



+ - 



dAp(A) log[sinh( 7 (A - i/2)) sinh( 7 (A + i/2))] 



(4.1) 



dAp(A)dA'p(A') log 



sinh( 7 (A - A' + i)) 
sinh( 7 (A- A')) 
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+ M- log sin 7 + log det D 
2 



The terms of order M 2 form the bulk free energy. In order to go further we have to 
analyze the behavior of the determinant. It is easy to see that the determinant of the 



matrix D : 
P(A), 

and therefore 



9<Pd 
9 V j 



is dominated by its diagonal elements; the latter are, by definition of 



dXj 2 



2tt— p(Xj) (4.2) 



M M f 

log det D « — log(TrM) H / dAp(A) log p( A) (4.3) 

2 2 J 

This gives us the expansion of the free energy up to linear terms in the size M. Note that 
similar expressions (with different densities p(A)) can be found for other states, as long as 
they have a proper M — > oo limit. 

As to the correlation functions, it is known that there is a general multiple integral 
representation for correlation functions in the thermodynamic limit ||17|| (see also [|18|, ||19| ). 
We shall not repeat the derivation here; let us simply mention that starting for example 
from (|3.4j ), one can prove the following formula |jT5] , pO|| : 

^ +1 ...^ + „)=2-»(^ n(n " 1)/2 rv.. d p„n . 5i h nl 7 (p, ~ pt > (4.4) 

V C / J-oo f* k smh 7(pj -pk-i) 



■j^j- sinh J 1 7(pj —i/2) sinh m 3 ^{pj + i/2) 



cosh m 7rp ? - 

This result is identical to the correlation functions of the 6-vertex model with periodic 
boundary conditions. 



Appendix A. A non-degeneracy property 

We assume in this appendix that q = e 17 is generic (i.e. not a root of unity, except the 
isotropic case q = —1). We also assume that the spectral parameters {p,k} do not form 
any "strings" (i.e. Im(fJ>k — Hi) 7^ 1 V/c, /). We consider two Bethe states |l) and |2) 
characterized by two sets {A], j = 1, . . . , m 1 } and {A^, j = 1, . . . , m 2 }. Bethe states are 
eigenstates of the set of commuting transfer matrices T(A), with corresponding eigenvalue 

T(A) \a) = 



Q a (X-i) 
Q Q (A) 



a(A) + 



Q a (A + z) 
Q a (A) 



d(A) 



(A.l) 
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where a(A) and d(A) are given by Eq. ( |2.13Q and are independent of the state, whereas Q a 
characterizes the {A"}: 

Q Q (A) = nsmh( 7 (A-A-)) (A.2) 

Note that the Bethe Ansatz Equations are simply the equations which ensure pole cancel- 
lation in the eigenvalue of T(A): Res T(A) |a), A _ AQ = 0. 

Because of the symmetry of the transfer matrices T(A) under the flip operator R, we 
are only considering states with m a < M/2. 

We now assume that |l) and |2^> have the same eigenvalue, that is 

~~ q^(a) + -QHAT d(A) = ~ww {) + ~wr ( } VA (A - 3) 

We rewrite this as 

a(A) [Q X (A + i)Q 2 (X) - Q 2 (A + OQ 1 ^)] = d(A) [Q 2 (A - i)Q 1 (A) - Q X (A - z)Q 2 (A)] 

(A ' 4) 

Up to an overall prefactor e - 2 ( M + m + m )i x ^ both left and right hand sides are polynomials 
in e 2lX of degree at most M + m 1 + m 2 . Furthermore they have the following 2M known 
zeroes: A = fik ± i/2, k = 1, . . . , M. If some coincide the zeroes have a multiplicity; 
however note that a //^ + i/2 cannot coincide with a \x\ — i/2 (since the are not allowed 
to form strings). There are now two situations: 

1) m 1 + m 2 < M. In this case we conclude directly that both sides of Eq. (A.4) are zero. 



2) m 1 + m 2 = M. Since m 1 < M/2 and m 2 < M/2, this can only happen if m 1 = 
m 2 = M/2. However in this case direct computation of the highest degree terms of the 
polynomials in Eq. ( [A .4 ) shows that they are zero, and therefore they are in fact of degree 



at most 2M — 1. Again this means that both sides of the equation are zero. 
In either case, we finally find 

Q^A + t) Q 2 (A + z) 

-QHAT = ~WW VA (A - 5) 

If q is not a root of unity, this implies immediately that {Aj} = {A|}. What we have 
proven is the following result: 

A.l. Two Bethe states with m < M/2 (i.e. S z > 0) correspond to the same eigenvalues 
of the T(A) (for all A) if and only if they are identical. 
This has the following two immediate corollaries: 
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A. 2. Two distinct Bethe states are orthogonal to each other. 

A. 3. A Bethe state with m = M/2 (i.e. S z = 0) is an eigenstate of the flip operator R. 

Remark: the situation is much more subtle if q is a root of unity. One way to see 
this is to consider Eq. ( |A.5|) , and assume now that q 2N = 1. One can add an extra "full 
string" of the form A; = Ao + ij, j = 1, . . . , N to one of the sets without modifying the 
corresponding Q function. This suggests extra degeneracy can appear when q is a root of 



unity between states with AS Z = N, which is precisely the phenomenon observed in [21 
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